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Abstract 

The extraction of meson-baryon scattering parameters from lattice QCD is complicated by the 
necessity, in some channels, of including annihilation diagrams. We consider a strategy to avoid the 
need for these extremely costly and noisy contributions. The strategy is based on simulations with 
an isospin chemical potential which, contrary to a baryon chemical potential, has no sign problem. 
When the isospin chemical potential is larger than a critical value, a charged pion condensate 
forms. Baryons propagating in the pion condensate will have their masses modified by pion- 
baryon scattering parameters. Consequently these parameters can be extracted from lattice QCD 
simulations in an isospin chemical potential, and we detail precisely which low-energy constants 
using baryon chiral perturbation theory. 
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I. INTRODUCTION 



In recent years, the determination of hadronic interactions via lattice QCD has shown 
substantial progress. In particular, calculations of meson-meson scattering lengths, such 
as that in the 1 = 2 channel of pion-pion scattering, agree to within a few percent of 
experiment 0, 0, !, H 9- Additionally within the past year, considerable progress has 
been made with multi- meson systems jsl, 0, 0] , two and three baryon interactions [1, T(J 11 



and several meson-baryon scattering processes 12| • The computation of a class of hadronic 
interactions, however, is hindered by the existence of annihilation diagrams, which, using 
current lattice methods, require all-to-all propagators. While some progress has been made in 
the computation of these propagators [131], the scaling to larger lattices is still prohibitively 
expensive for practical calculations of hadronic interactions. As a result, information on 
pion-nucleon scattering, and KN scatter ing, a p rocess thought to be important for kaon 



14 Il5l ; Il6l . 1 171 ll8| , has remained elusive to a first principles 



condensation in neutron stars 
lattice calculation. 

In this work, we propose a novel method for extracting the low-energy constants (LECs) 
for scattering processes with or without annihilation diagrams by employing an isospin 
chemical potential. In terms of lattice calculations, the addition of an isospin chemical 
potential results in a positive definite determinant, which avoids the fermion sign problem 
that plagues other finite density calculations, for more information and examples, see 



21, 22, 23 



An additional feature of an isospin chemical potential intrinsic to our analysis is 
the formation of a pion condensate when the chemical potential reaches a critical value [24], 
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This feature has been observed in previous lattice calculations [ 
the presence of this dynamically generated pion condensate, propagating baryons have their 
masses shifted by an amount depending on the value of the condensate, chemical potential, 
and a linear combination of LECs, with a dependence that can be calculated in heavy baryon 
chiral perturbation theory (HB^PT). As a result, precision measurement of baryon masses 
for several values of the chemical potential coupled with the HB^PT analysis gives an avenue 
for extracting scattering parameters, even for processes with annihilation diagrams. 

The organization of our work is as follows. In Sec. [TTJ the continuum QCD action with an 
isospin chemical potential is detailed along with the leading-order values for the condensates. 
In Sec. 1111} we present the two-flavor HBxPT analysis of nucleon and hyperon masses in the 
presence of the pion condensate. Our results show the linear combinations of LECs that 
can be extracted from doing several measurements at different chemical potentials above 
the critical point. In Sec. HVt the results of three-flavor HBxPT are derived by matching 
conditions to the two-flavor results. A discussion in Sec. |V] concludes our work. 



II. ISOSPIN CHEMICAL POTENTIAL 

Let us consider QCD with two light-quark flavors. The strange quark will be included 
below in Sect. IIVI Using a continuum Minkowski space action, the matter part of the QCD 
Lagrange density is 

c = i) 

where Jp is the <S77(3) C gauge covariant derivative, m q = diag(m, m) is the quark mass matrix 
in the isospin limit, and ip = (u, d) T an isodoublet of quarks. The isospin chemical potential, 
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Hi, appears as the time-component of a uniform gauge field. The term proportional to e is 
an explicit isospin breaking term. It is included here because the spontaneous breaking of 
isospin is an essential part of our analysis. Because symmetries do not break spontaneously 
at finite volumes, a small explicit breaking is necessary. This explicit breaking should go to 
zero as the size of the lattice goes to infinity. 1 

For vanishing quark masses, m — 0, and vanishing external sources, /// = e = 0, the QCD 
action has an SU{2) L ® SU{2) R symmetry that is spontaneously broken down to SU{2) V . 
For non-vanishing chemical potential, the QCD action has only a U(1)l <E> U(1)r symmetry 
associated with I 3 . This will be broken down to U(l)v by both spontaneous chiral symmetry 
breaking and the explicit symmetry breaking introduced by the quark mass. At finite \xi and 
finite e, the only continuous symmetry of C is the U(1)b associated with baryon number. 2 
We can write down an effective field theory of QCD that takes into account this pattern of 
spontaneous and explicit symmetry breaking. This effective theory is x?T, which is written 
in terms of the coset U G SU(2) L ® SU(2) R / SU(2) V) and at leading order has the form 

f 2 

£ = — [< D^UD^U 1 > +2A < M ] U + U ] M >] . (3) 
8 

Here the angled brackets denote flavor traces, and the external sources have been included 
in the terms 

M = s + ip, 

D^U = d^U + i[V^,U}, (4) 

3 

where the external vector potential is given by: = Hi^S^o, the pseusdoscalar source is 

2 

given by: p = e 2 ^-! an d the scalar source is just the quark mass, s = m q . 

The isospin chemical potential favors the presence of up quarks as opposed to down 
quarks. When it is larger than about m^, it becomes energetically favorable for the ground 
state to contain positive pions. Indeed, assuming a uniform value Uq for the vacuum ex- 
pectation value of U(x), we can minimize the effective potential to determine the vacuum 
(ground) state. Using the standard parametrization of SU(2), namely Uq = exp[ian ■ r], 
with n • n — 1, we find 

,. 21 h <m ". (5) 

exp[ia t J, \/jLi\ > m n 



1 In Euclidean space, the action can be written in the form Se = ^{'D + rn q )^p, where 

3 2 

V = ]p+ t i n J—+ie lh ^. (2) 

The Dirac operator T> satisfies the condition t 1 7 5 2?t 1 7s = T>\ which ensures that the product of u and 
d fcrmionic determinants is positive. 

With both and e non-vanishing there is no symmetry that prevents their renormalization. These 
parameters are only multiplicatively renormalized because, for example, the isospin chemical potential is 
the only term in the action that breaks C r J , and the explicit isospin breaking term is the only term in 
the action that breaks parity. Furthermore because we are interested in the limit of small e, the action will 
be approximately symmetric under the U(l)y associated with t 3 . Consequently m cannot be appreciably 
renormalized. 
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The angle a is determined, at lowest order in the chiral expansion, by the transcendental 
equation 

ml Ae 

cos a = — 5 ~ cot a, (6) 

Hi Hi 

where the pion mass satisfies the Gell-Mann-Oakes-Renner relation, m\ = 2Xm. For small 
values of e, one can determine a using a perturbative expansion of Eq. (jfJJ). This expansion, 
however, breaks down as one nears the critical value of the chemical potential from above. 
In the condensed phase, the condensates have values 

(■0V) = / 2 Acosa;, 
z( , 0t 2 75'0) = / 2 Asina. (7) 



III. BARYONS IN AN ISOSPIN CHEMICAL POTENTIAL: SU{2) 

When a hadron propagates in the background of a pion condensate its properties are 
changed due to the interactions of the hadron with the pions in the condensate. Their mass, 
for instance, is shifted by an amount closely related to the forward scattering amplitude of 
pions on the hadron. Unfortunately, the pions in the condensate are not exactly on-shell 
and no model independent relation can be found between pion-hadron scattering amplitudes 
and the mass shift. However, the chiral expansion of these two quantities are be related in 
the way described in this section. 



A. Nucleons 

We can address effects of the isospin chemical potential on nucleons by using the heavy 
nucleon chiral Lagrangian [32]. This effective theory accounts for spontaneous and explicit 
chiral symmetry breaking in the nucleon sector. The theory is written in terms of static 
nucleon fields, and thereby the nucleon mass can be treated commensurately with the chiral 
symmetry breaking scale. 

The leading-order heavy nucleon chiral Lagrangian is written in terms of the nucleon 
doublet, N = (p,n) T , and has the form 

£W = jVt (iv ■ D + 2g A S ■ A) N. (8) 

The vector is the nucleon four- velocity, while the vector, V^, and axial- vector, A^, fields 
of pions are defined by 

K = \ - ZrD%) . (9) 

The former appears in the chirally covariant derivative D^, which has the following action 
on the nucleon field 

(D^N), = d^Ni + ( WiVy (10) 

Appearing in the vector and axial-vector pion fields are £l, and £r which arise from the 
definition U = £l£r- Under a chiral transformation, (L,R) G SU(2) L <g> SU(2) R , the pion 
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fields transform as: — ¥ L ^l{x)U{x) , and £r(x) — > W {x)^r{x)W , where U{x) is the 
matrix entering the transformation of the nucleon field, Ni — > U(x)i : 'Nj. The left-handed 
and right-handed derivatives act according to the rules 

D*Zr = d^ R - i^. (11) 

For the case of an external vector field, the left- and right-handed vector sources coincide: 
L M = = Yfj,. In an isospin chemical potential, the spurions are given the final values 

U*) = Z(?)£o> (12) 

so that the vacuum value is Uq — £q. 

At leading-order using Lagrangian Eq. (jSJ), the shift in nucleon mass due to the chemical 
potential is the trivial result 

2 

where is the chiral limit value. Beyond leading-order, the nucleon mass receives cor- 
rections that depend on low-energy constants. These constants are the coefficients of terms 
in the second-order nucleon chiral Lagnrangian. Including relativistic corrections with fixed 
coefficients, the complete 3 second-order Lagrangian is 



M N = M$ - /x/cosa^-, (13) 



£(2) = jyt 



+ 2^F", s "] fin D u )-^{vA,S-D}+4 Cl <M + > 

N, (14) 



2M P '~ JL ^ m n 



+4 I c 2 - ) (v ■ A) 2 + 4c 3 A 2 + 4c 4 [£^ S V ]A„A V + 4c 5 M 



8i< 



where M+ = M+ - \ < M+ >, and M+ = \\(t} L M£} R + £ R M^ L ). The last term with 
coefficient C5 only contributes in the presence of strong isospin breaking. 

Utilizing the second-order nucleon Lagrangian, we arrive at the nucleon mass 



M 



N 



t 3 ( a 2 \ 

Mjy — hi cos a h 4ci (m 2 cos a + Ae sin a) + I c 2 f- + c 3 ) /jj sin 2 a (15) 

2 y 8M J 



The second-order correction is entirely isoscalar, and allows access to the low-energy constant 
Cj, and the combination c 2 + C3. One can vary the quark mass and isospin chemical potential 
to isolate these coefficients from the observed behavior of the nucleon mass in lattice QCD. 4 



3 When one considers external sources with non-vanishing field strength tensors, additional terms are 
present. The external sources we consider, however, are uniform in spacctimc. 

4 One can go further and impose isospin twisted boundary conditions on the quark fields, e.g. ip(x + Lz) = 
c~x.p[i0T 3 ]ip(x). This has the advantage of introducing non-vanishing spatial components of the vacuum 
value of A^, and leads to the ability to isolate more low-energy constants. For the example mentioned, 
the c 3 coefficient can be determined from the mixing angle between nucleons. 
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B. Hyperons 



Recently there has been interest in treating hyperons using SU (2) x?T [33|, |34j, |35|, 1361.137. 



38|. We consider hyperons in an isospin chemical potential, and determine which low-energy 
constants of SU(2) %PT can be determined from their masses. This is a natural approach 
because the isospin chemical potential produces only pion-hyperon interactions. 

In the stangeness 5 = 1 sector are the A and E baryons. The A baryon is an isosinglet, 
while the E baryons form an isotriplet. The leading-order Lagrangian in the S = 1 sector 
has the form 

£« = A+iu • <9A+ < E f (iv ■ D - A AE ) E > . (16) 
Here we have packaged the E fields into an adjoint matrix 

— E° E+ \ 

and use angled brackets to denote flavor traces. The covariant derivative acts as D^L = 
<9 M E + \Vft, E]. The parameter Aae is the SU{2) chiral limit value of the mass splitting, 
Aas — M^' — In writing the leading-order Lagrangian, we have omitted the axial 

couplings to pions as these terms are not relevant for our computation. Similarly we write 
down only the contributing terms from the second-order 5 = 1 Lagrangian 5 

L = At A [4cf < M+ > +24 < v ■ A 2 > +24 < A 2 >} 

+ < E f E > [4cf <M+> +24 < v ■ A 2 > +2cf < A 2 >] 

+4c 6 < E^ • A >< v ■ A E > +4c 7 < EU M >< A^Y, > . (18) 

Using the first- and second-order Lagrangians, we find the mass of the A baryon 

M A = M A 0) + 4cf {ml cos a + Ae sin a) + (c£ + <£)fi* sin 2 a, (19) 

and see that the low-energy constant , and the combination + can be accessed. For 
the E° baryon, we find the mass 

M E o = M§ ] + 4cf (ml cos a + Ae sin a) + (cf + cf sin 2 a. (20) 

At this order in the chiral expansion, the charged E baryons mix. We find the mass eigen- 
values 



M ± = M s o + (cg +C7)/z 2 sin 2 a±/i/ycos 2 a + (cf + cf ) 2 /i? sin 4 a. (21) 
The corresponding eigenstates are 

I") = sin||E + )-cos||E4, (22) 



To account for Kroll-Ruderman terms in the 5=1 
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r(0) 



M^ 0) ), and q 



|ffls/(^A T-^S ;j 1 '6 - ' L 6 12 

of the hyperon axial couplings and <7ae is that of [35 1 



sector. 

(o) 



replace — » 

f ^| s /(m| 0) 



c 2 



M^ -*). The normalization 
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where the mixing angle can be found from 



~ cos a .„„. 

cos 9 = ==. (23) 

a/ cos 2 a + (cf + C7 ) 2 /i 2 sin a 

Finally because the cf coefficients should scale with the inverse of the chiral symmetry 
breaking scale, the mixing between S baryons should be small for moderately sized 

In the strangeness S = 2 sector are the cascades. The case of the cascade is similar to the 
nucleon as the cascades transform as an isospin doublet, 5 = (H°,ET) T . Consequently the 
first- and second-order Lagrangians have the same form as those for the nucleon, Eqs. (|SJ) 

and ( fl4|) . One merely replaces the chiral limit mass of the nucleon, with the chiral 

limit mass of the cascade, M~ \ the nucleon axial charge, with the cascade axial charge, 
gss- Finally each of the low-energy constants Cj have different values for the cascades; these 
we denote by cj\ The cascade masses are thus given by 

M s = — cos a ~^+ 4c"f (ml cos a + Xe sin a) + ^cf ^foy + c t \ lA a - (24) 

Varying the quark mass and isospin chemical potential allows one to access cf, and the 
combination C2 + cf. 



IV. STRANGE QUARK AND SU(3) 

We can repeat the above analysis including a third light quark. The isospin chemical po- 
tential is added proportional to the A 3 Gell-Mann matrix, while the explicit isospin breaking 
term we choose in the direction of A 2 . For a massless strange quark, the QCD action with \ij 
and e both non- vanishing possesses a U(1)l <8> U(1)r symmetry associated with strangeness. 
The strange quark mass breaks this symmetry down to U(l)y- This symmetry precludes 
tilting the vacuum value of U in any strange direction. In SU(3), we must have 

Consequently the non-vanishing values of the condensates are 

~<^(1-V3A>> = A/ 2 , 

J(^(2 + v^A 8 )^) = A/ 2 cos a, 

i(^A 2 7 5 ^) = A/ 2 sin a. (26) 

It is now straightforward to determine the effects of the isospin chemical potential on the 
octet baryons. The octet baryons can be packaged into an SU(3) adjoint matrix B, given 
by 

7I S0 + 76 A S+ P \ 

B=\ S- -^o + ^A n . (27) 
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The relevant term of the leading-order 577(3) Lagrangian is 

£« =< B*[iV'D,B] >, (28) 

with D^B = d^B + [V M , B]. The kinetic term reproduces the isospin- dependent mass shifts 
for the octet baryons found above. 

The contributing terms of the second-order SU(3) chiral Lagrangian are 

£( 2 ) = 2b F < B ] [M + , B] > +2b D < B ] {M + , B} > +2b < B ] B >< M + > 

+h < B^, [A», B]] > +b 2 < B^[A^ {A*, B}] > +6 3 < B^{A^ {A*, B}} > 

+6 4 < B ] [v ■ A, [v ■ A, B]] > +65 < B ] [v ■ A, {v ■ A, B}] > 

+6 6 < B*{v ■ A,{v ■ A, B}} > +b 7 < B ] B >< A 2 > +b 8 < B ] B >< v ■ A 2 > . 

(29) 

Evaluating the masses of the octet baryons, we arrive at results identical to the SU (2) case 
but with SU(3) relations between parameters. The matching of parameters for the baryon 
masses have been given in [35j, while those for pion-baryon scattering have been given in [371 ] . 
Terms proportional to the strange quark mass give rise to shifts of the SU(3) chiral limit 
value of the baryon mass. Such terms are the leading contributions in SU(3) to the SU(2) 
chiral limit values of the baryon masses employed above. As our interest is with the /i/ 
dependence, we do not duplicate the matching relations between the chiral limit baryon 
masses here. Matching of Ai+ operators in the upper 2x2 block yields the relations 35 

4ci = b D + b F + 2b , 

4c^ = ^b D + 2b , 

4cf = 2b D + 2b , 

4cf = b D - b F + 2& , . (30) 

Matc hing of operators involving the product of two axial-vector pion fields gives the rela- 
tions [37] 

4(c 2 + c 3 ) = (h + 64) + {b 2 + 65) + (h + be) + 2{b 7 + b 8 ), 

4(4 + 4) = |(&3 + & 6 ) + 2(6 7 + 6 8 ), 

4(4 + cf) = 4(6 1 + 6 4 )+2(6 7 + 6 8 ), 
4(4 + 4) = -2(6i + b±) + 2(63 + b 6 ), 

4(4 + 4) = + h) - (62 + 6s) + (63 + h) + 2(6 7 + 6 8 ). (31) 

We have subsumed the Kroll-Ruderman terms here. From these matching conditions, we see 
that all three of the LECs for the quark mass dependent operators can be determined, along 
with four combinations of the eight double axial- insert ion operators. The overdetermination 
of these parameters, moreover, serves as a useful check on SU(3) symmetry. 



V. DISCUSSION 

Above, we show what aspects of meson-baryon scattering can be determined from sim- 
ulating lattice QCD with an isospin chemical potential. Our approach allows for the de- 
termination of the relevant LECs from first principles; and, specifically for pion-nucleon 
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processes, avoids the need to evaluate all-to-all propagators. Several attempts at extract- 
ing the non-derivative pion-nucleon couplings from analysis of experimental data have been 
attempted in the past. In f39[, nucleon-nucleon scattering data was used to determine the 
pion-nucleon constants through their effect on two-pion exchange contributions to the nu- 
clear force. Pion-nucleon scattering data has also been used for this purpose, for instance 
Significant uncertainties and discrepancies between different extractions remain. Our 



method can help resolve this situation. The values of the antikaon-nucleon couplings are 
considerably more uncertain, although attempts of extracting them from experimental data 



have been made, see, for example, [4ll . |42j |43| . |44j. Our framework can shed light on this 
issue by testing 577(3) symmetry predictions. 

Of particular interest is the coefficient ci, related to the nucleon a-term, which is im- 
portant for addressing dark matter interactions with nuclei, and has been the subject of 



numerous studies. The value estimated in |45| has been challenged [46]. The pion-nucleon 



constants also affect the dependence of the nucleon mass on the quark mass at fourth order 
in the chiral expansion, and this dependence has been used to extract them from lattice 
QCD calculations at different quark masses (for a recent extraction, see (13]). These ex- 
tractions are complicated, however, by the fact that the quark mass dependence of nucleon 
mass is barely compatible with chiral perturbation theory predictions at the current values 
of quark masses The method proposed here adds another tool that is somewhat 



complementary to the quark mass variation of the nucleon mass. In fact, in order to extract 
the pion-nucleon couplings one needs to know the nucleon mass for several values of the 
quark mass small enough so that chiral perturbation theory is trustable. But, apparently, 



this range might be small [49] . The dependence of the nucleon mass on the isospin chemical 
potential adds another handle that can be used for this purpose. The computational costs 
of changing the quark mass or the chemical potential are comparable. Of course, gauge 
configurations at different quark masses are also useful for other purposes besides nucleon 
mass extraction. The dependence on /i/, moreover, can also help in the extraction of other 
observables, making gauge configurations with fii more generally useful, and our proposal a 
realistic option. 
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